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SECTION  I 
INTRODUCTION 

Computational  aerodynamic  methods  that  are  significantly  better  than 
existing  techniques  are  needed  to  design  advanced  aircraft  configurations  in  a 
timely  and  cost-effective  manner.  The  need  for  Better  methods  is  due  in  part  to 
the  demand  for  computational  accuracy  brought  about  by  the  increased  aerodynamic 
sensitivity  at  transonic  speeds  of  advanced  technology  concepts  such  as  super¬ 
critical  wings,  active  controls,  variable  camber  wings,  and  laminar  flow  control. 
Also,  improved  aerodynamic  tools  are  required  to  handle  new  configurational 
concepts  such  as  (1)  high  aspect  ratio  wings,  spanloader  designs  with  thick  wings, 
and  winglets  for  transport-category  aircraft;  and,  (2)  swept  forward  wings,  var¬ 
iable  camber  wings  with  direct  lift  control,  canards,  and  blended-wing  concepts 
for  fighters.  Because  efficient  transonic  performance  continues  to  be  an  important 
design  requirement  for  military  and  commercial  aircraft,  there  is  also  a  particular 
need  for  improved  transonic  aerodynamic  computational  tools  to  accurately  calculate 
the  performance  for  each  of  these  new  technology  and  configurational  concepts. 

Significant  strides  have  been  made  in  the  development  of  3“D  transonic 
aerodynamic  design  and  analysis  codes  over  the  past  few  years.  Although  many 
of  the  methods  are  still  in  the  evolutionary  stage,  a  few  have  matured  to  the 
point  where  application  to  the  solution  of  aircraft  design  problems  is  practical. 

The  primary  objective  of  the  Advanced  Transonic  Technology  (ATT)  program 
was  to  demonstrate  that  performance  improvements  and/or  reduced  development 
time  and  costs  can  be  achieved  by  incorporating  3-D  transonic  methods  into 
aircraft  design  procedures.  The  attainment  of  the  objective  would  result  not 
only  in  very  efficient  advanced  technology  configurations  which  meet  future 
Air  Force  mission  requirements,  but  also  would  produce  an  experimentally 
verified,  efficient,  and  documented  transonic  configuration  design  method. 

The  overall  program  approach  and  the  individual  tasks  are  summarized  in 
Figure  1.  As  shown  in  that  figure,  the  program  was  conducted  by  performing 
six  technical  tasks  which  were  grouped  into  three  time  phases. 

During  Phase  I,  a  fighter  and  a  transport  configuration  employing  advanced 
technology  were  selected.  Also,  during  Phase  I,  the  3-D  transonic  analysis  codes 
were  evaluated,  and  necessary  modifications  for  design  procedure  use  were  made. 

Following  Air  Force  approval  of  the  selected  configurations,  computer  codes, 
and  the  design  procedure,  the  detailed  configuration  aerodynamic  designs  and 
wing  tunnel  tests  were  performed  in  Phase  II.  During  the  third  and  final  phase: 

(1)  wind  tunnel  data  were  used  to  determine  mission  performance;  (2)  theoretical 
predictions  and  experimental  results  were  compared;  and  (3)  a  final  production 
design  procedure  was  developed. 


SECTION  II 
DISCUSSION 

1.  TASK  1  -  CONFIGURATION  SELECTION 

During  Task  1,  a  transport  configuration  was  selected  as  the  baseline 
aircraft  for  the  detailed  aerodynamic  design  of  Task  3,  using  the  advanced 
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transonic  methods  developed  under  Task  2.  The  choice  of  the  baseline 
configuration  was  dictated  by  the  requirements  that  the  resulting  aircraft 
design  incorporate  advanced  technologies  and/or  unique  configuration  concepts. 
Additional  requirements  were  that  the  selected  configuration  offer  significant 
performance  improvements  for  an  existing  or  projected  Air  Force  mission,  A 
further  criterion  was  that  the  chosen  configuration  should  be  amenable  to 
design  using  the  new  3-D  transonic  methods, 

a.  Transport  Configuration 

There  were  two  possible  approaches  In  the  choice  of  the  transport 
configuration.  A  completely  new  design,  based  on  advanced  transport  studies 
being  conducted  at  Lockheed-Georgia,  could  be  selected;  alternatively,  the 
selected  configuration  could  be  based  on  a  refinement  of  an  existing  aircraft 
design.  The  political  and  economic  environment  relative  to  military  transport 
aircraft  at  the  time  of  the  initiation  of  the  contract  made  the  development 
of  a  modified  or  derivative  configuration  a  more  realistic  prospect  than  the 
development  of  a  completely  new  transport  aircraft,  where  design  requirements 
were  ill-defined.  Thus,  the  selected  transport  configuration  has  been  based 
on  the  redesign  of  an  existing  transport  configuration  with  increased  aero¬ 
dynamic  and  propulsive  efficiencies.  The  derivative  aircraft  exhibits  signi¬ 
ficantly  Improved  performance  for  Its  original  design  mission,  and  increased 
ability  to  perform  alternate  missions, 

Mission  Definition 

The  selected  configuration  is  a  derivative  of  the  C-141B  transport 
aircraft  designated  as  the  C-141B/AC2.  Figure  2a  shows  the  C-1A1B  config¬ 
uration  and  Figure  2b  shows  the  C-141B/AC2  configuration.  The  new  design  is 
intended  to  satisfy  long-range  cargo  transport,  tanker,  and  missile  launcher 
missions.  The  selected  configuration  is  designed  for  a  cruise  Mach  number  of 
0.80  with  a  payload  of  75,000  pounds  over  a  range  of  3,500  nautical  miles. 
Performance  constraints  include  a  field  length  of  7,500  feet,  and  an  initial 
cruise  altitude  of  35,000  feet. 

Configuration  Development 

The  primary  design  goal  in  the  development  of  the  C-141/AC2  design  is 
to  minimize  the  fuel  required  to  perform  the  design  mission.  The  C-141B/AC2 
configuration  was  sized  using  the  Lockheed-Georgia  General  Aircraft  Sizing 
Program  (GASP),  a  proprietary  computer  program  presently  used  for  all  company 
preliminary  design  studies.  The  program  accounts  for  the  interaction  of  the 
various  design  constraints  and  technical  disciplines  involved  in  the  aircraft 
design  process.  Drag  coefficients  and  weights  are  calculated  on  a  component 
basis  and  integrated  into  complete  aircraft  drag  and  weight.  The  propulsion 
system  is  sized  by  matching  cruise  thrust  requirements  within  the  constraints 
of  specified  takeoff  field  length. 

Technology  levels  for  the  various  disciplines  are  controlled  by  the  use 
of  input  adjusting  factors.  The  C— 1 4l B/AC2  aircraft  was  sized  using  advanced 
technology  levels  which  are  appropriate  to  an  Initial  operational  capability 
date  in  the  mid-1980's,  The  technologies  Incorporated  Into  the  design  are  as 
fol lows : 
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Aerodynamics  -  The  transonic  aerodynamic  technology  level  is  based  on  a 
method  of  Or.  Whitcomb  of  NASA-Langley  for  the  estimation  of  wing  thickness 
to  chord  ratio.  Thickness  is  calculated  as  a  function  of  lift  coefficient, 
cruise  Mach  number,  aspect  ratio  and  compressibility  drag  level  at  cruise. 
Airfoil  technology  level  is  selected  by  a  Mod  (drag  divergence  Mach  number) 
adjustment  which  has  been  assumed  to  be  zero  for  the  C-141B/AC2  1 985  IOC  date. 
Thus,  aerodynamic  estimates  used  in  GASP  for  this  study  are  conservative. 

Flight  Controls  -  The  C-141B/AC2  configuration  includes  automatic  load 
alleviation  of  wing  bending  loads,  relaxed  static  stability  and  a  stability 
augmentation  system.  The  wing  load  alleviation  system  employs  active  aileron 
and  elevator  controls  to  reduce  bending  loads  at  frequencies  up  to  the  first 
wing  bending  mode.  This  load  alleviation  permits  incorporation  of  a  higher 
aspect  ratio  wing  for  enhanced  fuel  efficiency.  The  relaxed  static  stability 
system  consists  of  automatic  flight  control  functions  which  provide  artificial 
static/dynamic  aircraft  stability  thereby  allowing  a  substantial  reduction  in 
tail  size.  Stability  augmentation  is  employed  in  the  lateral -d i rect ional  axes 
to  improve  flying  qualities. 

Structures  -  Composite  materials  are  used  for  the  secondary  structure, 
and  for  selected  reinforcement  of  the  primary  structure. 

Propulsion  -  The  C- 1 4 1 B/AC2  was  sized  using  current  technology  engine 
performance  represented  by  a  rubberized  version  of  the  CFM  56  engine. 

A  matrix  of  configuration  variables  was  input  to  the  sizing  program  to 
determine  the  minimum  fuel  aircraft  which  met  the  design  mission.  The  range 
of  variables  considered  was: 

Wing  loading  100  to  140  lb/ft2 

Aspect  ratio  8  to  12 

Wing  Sweep  10  to  25  deg 

Initial  cruise  altitude  31,000  tO  35,000  ft. 

Cruise  power  setting  0.7  to  1.0 

The  constraints  imposed  on  the  parametric  designs  were: 

Fuel  volume  ratio  >  1.0 

(Fuel  volume  available/ fuel  required  for  3500  NM) 

Field  length  <  7,500  ft. 

Aspect  ratio  <  12 

(To  avoid  excessive  wing  weights  or  the  necessity 
for  strut-braced  wings) 

Cruise  altitude  >  31,000  ft. 


3 


A  typical  sizing  plot  from  which  the  wing  planform  geometry  was  selected 
is  shown  in  Figure  3.  The  choice  of  wing  loading  and  sweep  is  established 
by  the  combined  requirements  of  minimum  mission  fuel,  and  the  fuel  volume  ratio 
required  to  meet  the  design  range  (FVR-1.0).  The  selected  wing  loading  of 
130  lb/ft2  and  wing  quarter  chord  sweep  of  25  degrees  satisfies  the  minimum 
mission  fuel  objective,  and  also  minimizes  the  aircraft  ramp  weight.  The 
selected  cruise  altitude  of  35,000  ft  and  a  wing  sweep  of  25  degrees  result 
from  the  parametric  variations  shown  in  Figure  A.  While  an  increase  in  initial 
cruise  altitude  above  35,000  ft.  would  save  more  fuel,  the  engine  thrust 
requirement  and  empty  weight  would  rise  significantly. 

Configuration  Characteristics  and  Performance 

The  resulting  C-1A1B/AC2  configuration  is  shown  in  the  general  arrangement 
drawing  of  Figure  2b.  A  complete  summary  of  the  aircraft  characteristics  as 
computed  by  the  GASP  Program  is  given  in  Figure  5.  The  major  features  of  the 
design  are: 

T.0.  gross  weight  261,100  lbs 

O.W.E.  118,300  lbs 

Block  fuel  61 ,500  lbs 

Cru i se  1 i ft  .58 

Start  cruise  wing  load  1 30  lb/sq  ft 

Wing  area  i960  sq  ft 

Quarter  chord  sweep  25° 

Aspect  ratio  12 

Taper  ratio  .4 

Average  thickness  10.9% 

Start  cruise  power  setting  80% 

Installed  thrust  38,800  lbs 

The  payload/range  capability  of  C- 1 4 1 B/AC2  aircraft  is  compared  to  the 
C- 1 4 1 A  and  C - 1 4 1 B  in  Figure  6.  The  C— 1 4 1 8/ AC2  OWE  is  about  80%  that  of  the 
C- 1 4 1 B ,  and  the  former  aircraft  requires  approximately  6A%  of  the  fuel  used  by 
the  latter  airplane.  Figure  7  shows  in  bar  chart  form,  the  impact  of  advanced 
technologies  on  the  C-141B/AC2  configuration  relative  to  the  baseline  air¬ 
craft  (C-141B). 

2.  TASK  2  -  SELECTION  AND  DEVELOPMENT  OF  METHODS 


The  objectives  of  this  task  were  to  evaluate  existing  three-dimensional 
transonic  methods,  improve  selected  methods,  and  develop  a  preliminary  version 
of  the  desired  transonic  design  method. 

The  improvement  of  3*0  transonic  codes  capable  of  treating  the 
sophisticated  configurations  developed  in  Task  1  was  critical  to  the  success 
of  the  design  method.  Numerous  3-D  codes  (e.g.  Refs.  1  through  5)  were 
available.  These  codes  differ  from  one  another  In  both  problem  formulations 
and  solution  techniques,  and  each  could  be  expected  to  excel  In  certain 


applications  and  to  be  less  successful  in  others.  Consequently,  more  than  one 
method  would  possibly  be  needed  in  the  design  procedure. 


We  emphasized  evaluation  of  the  NASA-Ames  3-D  transonic  small  disturbance 
code  because  this  code  is  more  easily  modified  to  provide  increased  configura¬ 
tion  flexibility  and  requires  less  computation  time  than  full  potential 
techniques. 

a.  Correlations  with  Transport  Test  Data 

Code  evaluation  for  transport  design  centered  on  correlations  of 
computations  using  the  Ames  3-D  transonic  code  with  C- 1 4 1  data  from  a  recent 
AEDC  16T  test  (Ref.  6).  A  Lockheed  version  (A7VP)  of  the  Ames  code  was  used 
in  those  calculations.  The  distinguishing  features  of  this  code  are  simplified 
input,  automatic  smooth  grid  generation,  a  revised  solution  loop  coding,  and 
control  stream  linking  with  boundary  layer  codes. 

The  first  solutions  were  computed  for  an  isolated  basic  C - 1 4 1  wing 
(pylon  and  nacelles  off)  using  theoretical  model  ordinates.  Inviscid  results 
computed  using  nonconservative  relaxation  (NCR)  and  fully  conservative  relax¬ 
ation  (FCR)  are  summarized  in  Figure  8.  The  following  observations  can  be 
made : 

1.  As  expected,  NCR  inviscid  results  are  in  better  (fairly  good) 
agreement  with  experiment  than  FCR  solutions. 

2.  The  gear  pods  (which  are  large  on  the  C- 1 4 1 )  have  a  significant 
influence  on  the  lower  surface  flow  as  evidenced  by  the  n  =  .19 
results.  This  disturbance  propagates  outboard  spreading  over  the 
ent i re  wing. 

3.  Leading  edge  pressures  are  not  well  predicted,  possibly  in 
part  due  to  inaccuracies  in  the  wind  tunnel  model  leading  edge 
contours. 

k.  Although  NCR  results  are  in  fair  agreement  with  experiment, 
inviscid  results  are  not  adequate  for  wing  design. 

5.  Fuselage  effects  (other  than  gear  pod)  are  not  of  critical 

importance  for  high  wing  transports  at  cruise  lift  coefficients, 
provided  the  fillet  is  well  designed. 

The  effects  of  viscosity  were  examined  using  the  Lockheed  control  stream 
linking  of  the  3-D  inviscid  transonic  code  with  a  strip  boundary  layer  routine. 

This  flexible  procedure  permits  the  use  of  any  boundary  layer  routine  and 
user  adjustment  of  fluid  wing  ordinates.  The  latter  feature  proves  to  be 
extremely  useful  in  modeling  the  trailing  edge  separation  predicted  on  the 
C- 1 41  wing. 

The  effects  of  viscosity  were  accounted  for  by  first  computing  NCR 
inviscid  potential  flow  solutions  because  such  solutions  are  a  better  approx¬ 
imation  to  the  real  viscous  flow  than  FCR  results,  and  use  of  NCR  inviscid 
pressures  in  the  subsequent  boundary  layer  analysis  should  yield  a  fairly 
good  approximation  of  the  viscous  flow.  Strip  boundary  layer  analyses  were 
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done  at  the  four  wing  defining  stations  (not  the  code  output  stations  shown 
in  Figure  8).  Trailing  edge  separation  was  predicted  at  each  control  station. 

Four  heuristic  models  were  used  in  the  separated  regions:  6*  »  constant, 
zf 1 u Id  airfoil  “  constant,  d6*/dx  -  constant,  and  dzfjuid  airfoil /dx  *  constant. 
The  last  model  produced  the  best  agreement  with  experiment.  The  modified 
5*'s  were  added  to  the  wing  ordinates  at  the  defining  stations  to  account  for 
boundary  layer  displacement  effects,  and  linear  lofting  was  used  between  the 
defining  stations  to  produce  the  "fluid"  wing.  The  inviscid  3-0  transonic 
code  was  then  used  in  the  FCR  mode  to  calculate  the  potential  flow  about  the 
"fluid"  wing. 

The  results  of  this  non-iterative  weak  interaction  viscous  analysis  are 
summarized  in  Figure  9.  Good  agreement  with  experiment,  especially  shock 
strength  and  location,  is  evident  in  these  results. 

The  results  shown  in  Figure  9  were  generated  using  measured  model  ordinates. 
Comparing  the  inviscid  NCR  solutions  with  those  of  Figure  8  (which  were  computed 
using  theoretical  ordinates)  show  little  difference  between  the  theoretical 
results  in  the  leading  edge  region  (where  differences  are  known  to  exist  between 
the  two  sets  of  ordinates),  and  generally  poor  agreement  with  experiment.  This 
disagreement  is  probably  a  shortcoming  of  small  disturbance  theory.  Conse¬ 
quently,  a  full  potential  equation  code  would  probably  be  needed  when  accurate 
leading  edge  solutions  are  required. 

b.  Code  Modifications 


Critical  code  modifications  to  permit  the  detailed  aerodynamic  design 
of  the  Task  1  configurations  were  the  addition  of  pylon/nacelle  modeling  for 
transport  design. 

The  ability  to  treat  pylon/nacelles  was  added  to  the  ATVP  3-D  transonic 
code.  This  was  done  in  two  phases.  In  the  first  phase,  the  objective  was  to 
produce  a  procedure  for  determining  the  effect  of  pylon/nacelle  geometry  on  the 
wing  flow  field,  and  not  for  computing  accurate  installed  pylon/nacelle  aero¬ 
dynamics.  In  this  case,  the  approach  selected  involves  first  the  calculation 
of  the  isolated  nacelle  flow  field.  This  was  done  using  a  modified  version 
of  a  full  potential  equation  code  written  by  Caughey  (Ref.  7)  for  axi symmetric 
nacel les. 

Once  the  nacelle  solution  is  calculated,  nacelle  boundary  conditions  are 
computed  on  the  surfaces  of  a  rectangular  prism.  The  prism  surfaces  correspond 
to  horizontal  and  vertical  grid  planes  in  the  3-0  transonic  code.  Nacelle 
boundary  conditions  are  imposed  on  these  grid  planes.  Small  disturbance  pylon 
boundary  conditions  are  incorporated  in  a  standard  manner  by  modifying  the 
difference  equations  to  include  the  pylon  surface  slopes,  dypylon^x»  'n  tbe 
spanwise  difference  terms. 


In  the  second  phase,  an  extended  small  disturbance  flow  through  nacelle 
code  was  written  and  linked  to  the  ATWP  3-D  wing  code.  This  combined  code  (TALA) 
can  treat  mildly  non-ax 1  symmetric  nacelles  at  small  angles  of  attack  and  yaw.  It 
also  cycles  back  and  forth  between  nacelle  solutions  and  wing  solutions  to  obtain 
the  mutual  Interference  effects  between  nacelle  and  wing. 
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c.  Design  Procedure 

The  transport  wing  design  procedure  is  shown  schematically  in  Figure  10. 
The  procedure  is  based  upon  the  use  of  an  isolated  wing  code  for  the  transonic 
design.  Interference  pressures  are  obtained  from  a  transonic  wing-pylon-nacelle 
code,  and  also  from  a  more  economical  subsonic  panel  method  (for  lower  surface 
gear  pod  interference). 

The  key  element  in  the  procedure  is  the  use  of  numerical  optimization  in 
the  wing  design.  The  wing  design  code  was  developed  by  linking  Vanderplaats 
constrained  function  minimization  program  (CONMIN,  Ref.  8)  with  the  ATWP  3-0 
wing  code.  As  part  of  this  task,  a  design  objective  and  appropriate  constraints 
were  selected,  and  the  permissible  geometric  perturbations  (i.e.  design  vari¬ 
ables)  were  specified.  The  entire  process  is  described  in  detail  in  the 
following  paragraphs. 

Design  Objective  and  Constraints 

To  avoid  using  inaccurately  calculated  quantities  such  as  drag  in  the 
optimization  procedure,  the  program  was  developed  to  permit  the  design  of  wings 
with  specified  chordwise  pressure  distributions.  The  capability  of  examining 
two  pressure  design  objectives  was  provided.  One  design  objective  was  the 
minimization  of  the  RMS  deviation  between  the  target  and  actual  pressures: 


OBJ,-  [l(Cp  -  CpD)2/NJ  i 

where  N  is  the  number  of  pressure  coefficients,  and  CpD  is  the  target 
pressure  coefficient.  The  second  objective  considered  was 


0BJ2  -  (Cp  -  Cp0) ;  Constraint  -  CpD  <  Cp* 

Notice  the  constraint  required  to  make  the  second  objective  meaningful. 

Both  objectives  were  tried  in  the  design  study,  and  the  first  objective 
proved  to  be  superior.  Consequently,  OBJ,  was  used  in  the  transport  wing 
design. 

Design  Variables 

Consistent  with  established  wing  geometry  definition  procedures,  the 
wing  geometry  is  determined  by  specifying  the  airfoil  sections  at  various 
geometric  control  span  stations  and  connecting  these  sections  by  linear  loft 
elements.  At  each  control  station,  the  permissible  surface  perturbations 
are  listed  in  Table  1.  The  magnitude  of  each  of  these  1 A  surface  perturbation 
functions  plus  the  section  twist  angle  are  the  fifteen  design  variables  avail¬ 
able  for  each  surface  at  each  geometric  control  span  station.  Thus,  for  a 
four  control  station  wing,  if  all  the  surface  perturbations  were  used,  and 
if  all  the  sections  were  designed  simultaneously,  a  total  of  15  variables  per 
surface  per  station  x  2  surfaces  x  k  stations  -  120  design  variables,  would 
be  required. 


7 


Implementation 

The  simultaneous  use  of  120  design  variables  would  result  in  an 
inordinately  long  computer  run  (greater  than  10  hours  on  a  CDC  7600).  Job 
turn-around  time  on  such  a  run  would  be  very  long,  and  a  mistake  would  waste 
a  lot  of  computer  time.  Consequent! y ,  wing  design  was  accomplished  in  a  series 
of  steps.  First,  the  upper  surface  was  designed  one  section  at  a  time  pro¬ 
ceeding  from  the  root  to  tip.  Next,  the  lower  surface  was  designed,  again 
proceeding  from  root  to  tip.  The  optimization  is  done  using  ATWP  in  the  fully 
conservative  relaxation  (FCR)  mode,  and  the  desired  viscous  pressure  distribution 
is  specified.  Consequently,  the  design  procedure  produces  the  "fluid"  wing 
geometry  (that  is,  the  desired  solid  wing  geometry  plus  the  boundary  layer 
displacement  thickness).  The  fluid  wing  is  then  analyzed  using  FCR  and  the 
entire  process,  or  parts  thereof,  are  repeated  as  required  to  produce  the 
desired  pressures. 

Extraction  of  Solid  Wing  Geometry 

The  wing  contours  produced  by  the  optimization  include  the  boundary  layer 
displacement  thickness,  6*.  The  solid  wing  geometry  is  found  by  subtracting 
6*  from  the  computed  wing  contours  at  each  of  the  design  stations. 


TABLE  1  TRANSPORT  WING  DESIGN  VARIABLES 

.25  -20  (X/C) 

V(1)  -  3.89  (X/C)  (1  -  X/C)e 

V (2 )  -  10.68  (X/C)'  (1  -  X/C)e  U/C; 

V(j)  -  sin3  (ir(X/C)rj),  j  -  3,  12 

V(13)-  (X/C)8 

V(1A)  -  (X/C)20  ea  /  1  -  (X/C) 

a  -  .5/(1  -  (X/C)M)  -  20/ (X/C) ^ 
(X/C)M  h(X/C)  for  max  camber 
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TABLE  1  TRANSPORT  WING  DESIGN  VARIABLES  (continued) 

Sine  deformations 


J. 

r  j 

<X/C»H,x  , 

3 

.231 

.05 

4 

.301 

.10 

5 

.431 

.30 

6 

.576 

.30 

7 

.756 

.40 

8 

1.000 

.50 

9 

1.357 

.60 

10 

1  -943 

.70 

11 

3.106 

.80 

12 

6.579 

.90 

A  conventional  (Ref.  9)  2-D  integral  boundary  layer  code  is  used  to  compute 

6*. 

Analysis  of  Optimized  Wing 

The  performance  of  the  optimized  solid  wing  is  then  investigated  using  both 
the  full  potential  and  extended  small  disturbance  transonic  codes.  This  analysis 
should  use  the  viscous  linking  technique  of  the  design  method,  and  a  coupled  viscous 
analysis  technique. 

Second  Iteration  Wing  Optimization 

If  analyses  show  that  the  wing  does  not  perform  as  expected,  a  second 
pass  is  made  through  the  design  process.  All  or  part  of  the  procedure  may 
be  repeated  as  dictated  by  the  cause  of  the  unacceptable  performance. 

3.  TASK  3  -  CONFIGURATION  DESIGN  AND  WIND  TUNNEL  TEST  -  TRANSPORT 

The  design  procedure  developed  in  Task  2  was  used  to  design  a  new  wing 
for  the  C- 1 4 1 B/AC2  configuration  selected  during  Phase  I.  The  goal  of  this 
study  was  to  improve  the  wing  aerodynamics  and  at  the  same  time  increase  the 
wing  thickness  for  the  Mach*. 30,  Cl  *.60  design  condition.  Increased  wing 
thickness  was  sought  to  increase  the  fuel  volume,  and  to  reduce  wing  weight. 

The  new  wing  was  built  and  tested  in  the  Lockheed-Georg la  Company's  compress¬ 
ible  flow  transonic  wind  tunnel  (CFWT). 
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a. 


Starting  Wing  Selection 


Two  starting  wings  were  considered.  One  wing  used  airfoil  sections 
developed  at  NASA-Langley  as  part  of  the  Energy  Efficient  Transport  Program 
(Ref.  10).  The  second  wing  used  airfoils  designed  at  Lockheed -Geo rg i a  in 
support  of  in-house  configuration  Advanced  Transonic  Aircraft  (ATA)  design 
studies  (Ref.  11).  Two  wings  were  analyzed  using  ATWP;  results  are  summarized 
in  Figure  11.  These  results  show  that  both  wings  have  well-behaved  design 
point  pressures.  Thus,  either  wing  is  satisfactory  for  starting  the  design 
process. 

The  wing  with  the  Lockheed  airfoils  was  selected  as  the  starting  wing 
because  the  airfoils  were  systematically  developed  using  state-of-the-art 
airfoil  inverse  (Ref.  12)  and  analysis  (Ref.  13)  transonic  methods.  Consequently, 
this  initial  wing  design  lends  itself  well  for  inclusion  in  a  formal  design 
procedure. 

b.  Interference  Pressures 

Using  in  isolated  wing  code  during  numerical  optimization  is  a  key  feature 
of  the  design  method  because  this  minimizes  computer  resource  requirements. 

During  Ph«i*  I  the  isolated  wing  version  of  ATWP  was  shown  to  yield  good  pre¬ 
dictions  oi  uppe>  surface  pressures  for  high-wing  configurations.  The  inter¬ 
ference  pr-ssure  perturbations  are  confined  to  the  lower  surface  and  are  pri¬ 
marily  du^  o  *he  gear  pod.  Since  the  flow  is  subcritical  on  the  lower  surface, 
and  S'vce  the  interference  pressures  are  small,  they  can  be  calculated  using 
any  subsonic  panel  method  which  provides  very  good  geometric  resolution  in  the 
wing-body-gear  pod  area  as  shown  in  Figure  12. 

The  calculated  isolated  wing  and  wing-body-gear  pod  pressures  are  compared 
in  Figure  13  where  the  lower  surface  interference  pressures  are  apparent.  These 
interference  pressures  are  taken  into  account  in  the  specification  of  the  target 
Isolated  wing  design  pressures.  The  interference  effects  of  the  pylon-nacelle 
were  determined  by  including  pylon-nacelle  analysis  capability  in  the  ATWP  3-D 
transonic  analysis  code. 

c.  Design  Pressures 

Target  wing  pressures  were  originally  specified  near  the  wing  root,  break 
and  tip.  The  upper  surface  pressures  were  selected  to  provide  a  weak  shock 
wave  on  the  outer  panel  near  the  S0%  chord  station.  The  root  pressures  were 
selected  to  minimize  isobar  unsweeping  and  to  avoid  large  trailing  edge  pressure 
gradients  which  might  result  in  the  formation  of  a  strong  trailing  edge  shock 
wave. 


d.  Numerical  Optimization 

The  numerical  optimization  procedure  described  in  Subsection  2-C  was  used 
to  determine  the  wing  geometry  which  produces  the  desired  wing  pressures.  The 
wing  geometry  was  determined  by  designed  the  four  wing  control  stations  shown 
in  Figure  lA  and  using  linear  lofting  to  generate  intermediate  ordinates. 

Early  numerical  experiments  showed  that  the  additional  n  ■  .6  design  station 
was  needed  because  perturbations  of  the  wing  tip  section  were  inadequate  to 
control  mid-semispan  and  tip  pressures  simultaneously.  A  constant  normalized 
section  wing  carry-through  was  used. 
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Of  note  also  in  Figure  14  is  the  specification  of  tip  pressures  near  the 
n  -  .85  station.  This  choice  was  made  because  of  the  inaccuracy  of  computed 
results  near  the  wing  tip. 

Extended  Small  Disturbance  Design  -  The  ATWP  extended  small  disturbance 
program  was  used  in  the  initial  wing  design.  The  final  wing  pressures  are  com¬ 
pared  with  the  target  pressures  in  Figure  15.  Also  shown  here  is  the  agreement 
between  target  and  computed  pressures  (labeled  as  upper)  after  the  upper  surface 
design  of  each  span  station.  The  agreement  between  target  and  computed  pressures 
is  good.  The  computed  pressures  do  produce  a  slightly  stronger  than  desired 
shock  wave. 

The  solid  wing  geometry  was  then  analyzed  using  both  full  potential  (FL022, 
Ref.  14)  and  extended  small  disturbance  (ATWP)  codes  with  weak  interaction 
viscous  linking.  The  results  of  these  calculations  are  shown  in  Figure  16.  On 
the  premise  that  the  FPE  results  are  correct,  these  data  show  that  the  ESO  code 
entirely  mispredicts  the  wing  leading  edge  flow  field,  and  this  error  causes 
complete  disagreement  between  ESD  and  FPE  results. 

Full  Potential  Design  -  The  failure  of  the  ESO  optimization  to  accurately 
design  the  wing  leading  edge  made  a  second  pass  through  the  design  procedure 
using  a  full  potential  equation  analysis  code  necessary.  The  optimization  code 
was  modified  for  this  purpose,  and  the  FPE  optimization  was  done  using  the  same 
design  variables  and  objectives  used  in  the  ESO  optimization.  However,  the 
target  pressures  for  the  FPE  design  were  modified  to  produce  a  slightly  weaker 
shock  wave. 

The  span  load  distribution  of  the  resulting  wing  was  not  satisfactory 
because  the  outboard  section  of  the  wing  was  too  highly  loaded,  while  the 
root  section  was  unloaded.  The  lower  surface  of  the  root  was  modified  using  a 
simple  Lockheed  linear  design  method  to  increase  the  section  loading  at  a  small 
sacrifice  in  wing  root  thickness.  The  wing  twist  was  also  adjusted  using  a 
very  efficient  panel  method  to  decrease  the  tip  loading. 

The  resulting  wing  was  analyzed  using  both  the  FPE  and  ESD  codes  with 
viscous  linking  (as  described  in  Subsection  2-a) .  The  results  are  summarized 
in  Figure  17.  The  codes  produce  results  in  good  agreement  with  one  another, 
indicating  that  ESD  methods  yield  accurate  results  if  the  leading  edge  is 
properly  designed.  The  wing  pressures  are  quite  satisfactory.  There  is  no 
tendency  for  isobars  to  coalesce  near  the  root  trailing  edge,  nor  is  there  a 
tendency  for  a  leading  edge  shock  wave  to  form.  The  desired  mid-chord  shock 
is  weak  (normal  Mach  number  less  than  1.16).  Consequently,  this  FPE-designed 
wing  was  selected  as  the  final  C-141B/AC2  wing  design. 

e.  Wind  Tunnel  Test 


Models 

An  existing  .0188  scale  C-I4IA  semi-span  model,  was  used  to  obtain 
baseline  data  for  comparison  with  the  C-141B/AC2  configuration.  The  wing 
on  this  model  Is  instrumented  with  126  surface  static  pressure  taps  located 
at  three  span  stations.  The  new,  .0188  scale  C-141B/AC2  wing,  was  machined 
from  a  solid  billet  of  17-4PH  stainless  steel  by  a  numerically  controlled 
milling  machine.  The  model  was  hand  finished  to  a  tolerance  of  ±  .002  in. 
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Actual  model  coordinates  were  measured  at  several  stations  by  a  precision 
measuring  machine,  and  compared  to  lofted  coordinates  used  to  define  the 
wing.  The  C-141B/AC2  wing  had  a  total  of  1 40  static  pressure  taps  located 
in  chcrdwise  rows  at  5  span  stations.  The  upper  surface  pressure  orifices 
were  installed  by  drilling  completely  through  the  wing  so  that  all  pressure 
tube  routing  is  on  the  wing  lower  surface.  This  installation  technique 
makes  the  upper  surface  completely  free  of  tube  routing  channels  which  can 
cause  possible  discontinuities  in  curvature  that  may  affect  the  supersonic 
flow  on  the  upper  surface. 

The  fuselage  used  in  these  tests  was  a  C- 1 4 1 A  model  fuselage  originally 
built  for  the  C-141  semi-span  wing.  Using  the  A  model  fuselage  would  make 
little  difference  aerodynamical ly  and  would  save  the  cost  of  fabricating  a 
new  B  model  fuselage.  A  sketch  of  the  fuselage  is  shown  in  Fig.  18.  The 
C-141  wing  planform  is  shown  in  Fig.  19*  Table  2  lists  the  C-141  wing 
pressure  tap  locations.  Fig.  20  shows  the  C-141B/AC2  wing  planform.  Table 
3  lists  the  AC2  pressure  tap  locations.  Table  4  lists  the  AC2  wing  ordinates 
as  measured  by  a  numerical  measuring  machine.  Fig.  21  shows  a  direct  compar¬ 
ison  of  the  C-141  and  C-141B/AC2  wing  planforms. 

Teat  Facility 

The  general  arrangement  of  the  Lockheed  Compressible  Flow  Wind  Tunnel 
(CFWT)  is  shown  in  Fig.  22.  The  tunnel  is  of  the  blow-down  type,  exhausting 
directly  to  the  atmosphere.  The  air  storage  capability  is  13,000  ft3  at 
600  PSIA.  A  sleeve  type  control  valve  accurately  maintains  the  settling 
chamber  stagnation  pressure  at  selected  pressure  less  than  or  equal  to  the 
250  psia  maximum  at  a  mass  flow  rate  less  than  2400  Ib/sec.  The  test 
section  is  20.0  In.  wide  by  28.0  in.  high  by  72.0  in.  long  and  is  enclosed 
in  a  12.0  ft.  diameter  plenum  chamber.  The  top  and  side  walls  of  the 
three-dimensional  test  section  have  variable  porosity  capability  (from  0  to 
10  percent),  obtained  by  sliding  two  parallel  plates  with  .250  in.  diameter 
holes  slanted  60  degrees  from  the  vertical.  The  bottom  wall,  where  the 
model  is  mounted,  is  not  porous.  The  model  is  mounted  on  a  five-component 
balance  located  in  the  floor.  The  balance  and  model  rotate  together  on  a 
turntable  to  vary  angle  of  attack.  A  bleed  duct  is  located  21  in.  ahead 
of  the  balance  centerline  to  remove  the  wind  tunnel  boundary  layer.  The 
boundary  layer  bleed  system  has  an  independent  control  valve  and  exhausts  to 
atmosphere  through  a  separate  pipe  system.  The  main  features  of  the  test 
apparatus  are  illustrated  in  Fig.  23.  A  more  detailed  description  of  the 
facility  may  be  found  in  Reference  15. 

Ine trumen taticm 

Instrumentation  for  this  test  program  consists  of  a  five-component 
strain  gauge  balance  to  measure  model  aerodynamic  forces.  Six  far-field 
pressure  rails,  each  containing  31  static  pressure  taps,  were  mounted  on 
the  tunnel  walls  as  shown  in  Fig.  23  and  detailed  in  Fig.  24.  The  far-field 
measurements  were  extended  to  the  symmetry  plane  by  a  row  of  fourteen  stream- 
wise  pressure  taps  along  the  tunnel  floor  on  each  side  of  the  model.  These 
pressure  orifices  were  displaced  2.0  in.  inward  from  the  tunnel  sidewalls. 


Measurements  of  the  static  pressures  on  the  wing  surface  and  the  wall 
rails  were  made  using  electronical ly  actuated  pressure  scanning  valves. 

The  full-scale  range  of  the  quarter  percent  accuracy  Statham  pressure  trans¬ 
ducers  in  the  scanning  valves  were  selected  to  provide  maximum  accuracy  for 
the  wind  tunnel  conditions  tested  (wail  rails  ±12.5  psl  and  airfoil  pressures 
±50  psi).  CEC  force  balance  pressure  transducers  were  used  in  conjunction  with 
CEC  servo  amplifiers  to  provide  a  precise  measurement  of  the  atmospheric 
pressure,  stagnation  pressure,  and  test  section  static  pressure  to  0.05%  of 
the  250  psi  capacity.  These  transducers  allow  determination  of  the  test 
section  Mach  number  to  an  accuracy  of  ±.002  at  the  highest  stagnation  pressure. 

Angle  of  attack  was  measured  with  a  calibrated  potentiometer  operated 
by  the  angle-of-attack  drive  mechanism. 

The  balance  used  to  measure  forces  on  the  model  was  designed  and  built 
by  Lockheed-Cal i forn ia  Co.  It  is  a  strain  gauge  five -component  temperature 
compensated  balance  designed  to  the  following  load  capacities: 


Normal  force 
Axial  force 
Pitching  moment 
Rol 1 ing  moment 
Yawing  moment 


±750  lbs. 

±75  lbs. 

±1800  in.  lb. 
±6000  in.  lb. 
±600  in.  lb. 


Design  accuracy  of  the  balance  is  ]/k%  of  design  load  and  \/2%  of  applied 
load.  Raw  pressure  and  balance  data  were  recorded  on  magnetic  tape  utilizing 
the  CFWT  high  speed  data  acquisition  system.  The  data  acquisition  system 
consists  of  a  Lockheed  Electronics  Compay  MAC-16  computer  and  associated 
peripheral  equipment.  The  raw  data  was  reduced  to  coefficient  form  with  a 
TI-990  computer.  Machine  plots  were  made  on  a  Calcomp  765  plotter. 


Data  Reduction 


Balance  data  -  The  force  balance  data  were  acquired  by  taking  one 
hundred  samples  over  a  one  second  time  span.  The  average  of  the  100  samples 
was  reduced  to  pounds  of  force  and  inch  pounds  of  moment  by  multiplying  the 
five  measured  values  by  the  calibration  matrix.  The  resulting  measurements 
were  converted  to  conventional  coefficient  form  using  the  tunnel  conditions 
at  the  time  the  force  measurements  were  made. 


Pressure  data  -  All  static  pressure  measurements  from  the  wing  and 
far-field  rails  were  reduced  to  standard  coefficient  form  using  instan¬ 
taneous  tunnel  conditions  at  the  time  each  pressure  was  measured.  Pressure 
coefficients  for  each  spanwise  station  were  numerically  integrated  to  determine 
a  local  lift  coefficient. 

(2)  Tests  and  Methods 
Test  Conditions 


r  atm qttematoiWfcW.Mi.  ■ 
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The  configurations  were  tested  over  a  wide  range  of  conditions  to  pro¬ 
vide  data  at  off-design  as  well  as  design  conditions.  The  test  Reynolds 
number  based  on  mean  aerodynamic  chord  was  nominally  5  million.  Although 
the  wind  tunnel  was  capable  of  much  higher  unit  Reynolds  numbers,  the  actual 
test  values  were  limited  by  the  capacity  of  the  wind  tunnel  force-balance. 

The  angle  of  attack  and  tunnel  Mach  number  were  varied  over  a  wide 
range:  from  zero-lift  to  above  design  values  and  subcritical  speeds  through 
arag-rise,  respectively.  A  complete  summary  of  the  test  matrix  is  given  in 
Table  5.  Table  6  gives  a  detailed  breakdown  of  the  run  numbers  corresponding 
to  each  test  condition  for  the  C— 1  ^4 1 B/AC2  wing.  Table  7  gives  the  run  numbers 
for  the  C- 1 4 1 B  wing. 

Most  of  the  wind  tunnel  testing  was  conducted  at  a  fixed  wall  porosity 
of  because  previous  tests  in  the  CFWT  had  indicated  this  value  would 
provide  minimum  wall  interference  effects. 

Transition 

For  the  unit  Reynolds  numbers  of  these  tests,  experience  with  the  CFWT 
facility  indicates  that  significant  regions  of  laminar  flow  would  exist  on 
a  smooth  model.  Therefore,  to  fix  transition  and  to  eliminate  the  uncer¬ 
tainties  of  the  transition  location,  a  transition  strip  was  applied  to  each 
model.  The  transition  strips  were  applied  to  both  the  upper  and  lower 
surface  and  were  located  and  sized  according  to  the  guidelines  of  reference  16. 
The  strips  were  thus  located  a  fixed  distance  from  the  leading  edge  equal  to 
5%  of  the  mean  aerodynamic  chord.  The  grit  consisted  of  .0023  in.  diameter 
glass  beads,  set  in  a  lacquer  fixative.  The  width  of  the  strips  was  constant 
at  .05  in.  Strips  were  also  placed  on  the  fuselage. 


Boundary  Layer  Removal  System 

One  unique  aspect  of  the  test  facility  used  in  this  test  program  was  the 
bleed  system  used  to  remove  the  wind  tunnel  boundary  layer  ahead  of  the  model 
(Fig.  23).  The  operational  procedure  used  in  all  of  the  testing  was  to  set 
the  bleed  rate  so  that  the  static  pressures  measured  at  the  leading  edge  of  the 
flat  plate  which  formed  the  top  of  the  bleed  duct  did  not  indicate  any  angle- 
of-attack  loading. 

Tunnel  Calibration 

Fig.  25  shows  that  the  CFWT  Tunnei  Mach  number  distribution  is  uniform. 
The  force  balance  used  for  these  tests  was  calibrated  prior  to  the  conduct 
of  the  tests. 

Fuselage /Wing  Seal 

The  wing/fuselage  configurations  were  tested  with  the  fuselage  non¬ 
metric.  The  fuselage  was  cut  out  so  that  the  wing  could  fit  through  the 
fuselage  with  a  nominal  gap  of  .040  in.  This  gap  was  filled  with  a  foam 
seal.  Any  possible  tare  effects  of  the  foam  seal  on  the  balance  reading 


were  accounted  for  by  applying  loads  to  the  clean  wing,  and  then  to  the  wing  + 
seal  +  fuselage  to  obtain  corrections  to  the  balance  readings.  In  any  case 
required  corrections  turned  out  to  be  small. 

1*0  Test  Results 

We  will  show  plots  of  all  of  the  force  data  obtained,  and  selected 
samples  of  the  pressure  data.  All  of  the  data  from  this  test  can  be  obtained 
on  magnetic  tape. 

Force  Data 

AH  of  the  force  data  obtained  during  the  wind  tunnel  test  is  shown  in 
igures  26  through  32.  The  data  is  grouped  by  configuration,  and  three 
figures  are  shown  for  each  configuration  (CL  vs.  a,  CL  vs.  CM,  and  CL  vs.  CD) . 


These  force  data  plots  will  be  discussed  in  Subsection  4.  jhe  force 
data  for  the  baseline  C— 1 4 1 B  wi ng-body-gearpod  configuration  and  the  C -*  1 4 1 B/ 

AC2  wi ng-body-gearpod  configuration  were  used  to  determine  performance  im¬ 
provements.  These  performance  comparisons  will  be  discussed  in  detail  in  section 
Subsection  5. 

Pressure  Data 

Much  of  the  pressure  data  obtained  for  the  C-141B/AC2  wing  along  con¬ 
figuration  is  given  in  the  discussion  of  code  correlations  later  in  this 
report.  Here,  we  will  show  only  some  sample  data. 

Fig.  33  is  presented  to  show  the  effects  of  the  gearpod,  on  lower 
surface  pressures.  These  effects  are  apparent,  especially  near  the  inward 
part  of  the  wing. 

Fig.  34  is  a  set  of  isometric  plots  of  the  pressure  data  for  the  C-141B/ 
AC2  wing  alone  configuration  at  Mach  0.8,  At  the  lower  angles  of  attack,  the 
outboard  wing  section  shows  negative  lift  in  the  first  20  to  30  percent  chord. 
The  shock  wave  at  approximately  60*  chord  is  apparent.  This  shock  wave  was 
specified  in  the  design. 
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4.  TASK  4  -  AERODATA  EVALUATIONS 


The  purpose  of  Task  IV  was  to  evaluate  the  wind  tunnel  data  and  compare 
this  data  with  computed  results. 

a.  Force  Data 

All  of  the  force  data  obtained  during  the  wind  tunnel  test  is  shown 
in  Figures  26  through  32.  In  this  section  we  will  discuss  that  part  of 
the  force  data  which  is  relevant  to  a)  configuration  performance  comparisons, 
and  b)  code  correlations. 

Figures  27a  to  27c  are  plots  of  vs.  a,  vs.  CM,  amd  C(_  vs.  Cp, 
respectively,  for  the  C- 1 4 1 B/AC2  wing  alone  configuration.  For  the  Cl  vs.  o 
curves  in  Figure  27a,  as  Mach  number  increases,  the  slope  becomes  steeper 
in  the  mid  Cl  range.  At  higher  values  of  lift  and  Mach  number,  the  slope 
decreases,  and  the  data  become  non-linear,  reflecting  the  effect  of  super¬ 
critical  flow  on  the  wing.  The  pitching  moment  curves  in  Figure  27b  are 
generally  non-linear.  Additionally,  the  pitching  moments  are  generally 
positive,  whereas,  for  most  transport  type  wings,  the  pitching  moments 
would  be  expected  to  be  negative.  This  can  be  explained  by  the  fact  that 
for  this  wing,  the  break  occurs  on  the  leading  instead  of  the  trailing  edge. 
As  a  result,  the  inboard  part  of  the  wing  makes  a  large  positive  contribution 
to  the  total  moment  about  the  quarter  chord  point  of  the  mean  aerodynamic 
chord  (MAC).  The  drag  polars  (Figure  27c)  vary  as  expected,  with  drag  in¬ 
creasing  for  a  fixed  lift  coefficient,  as  the  Mach  number  increases. 

Figures  29a  through  29c  give  the  Cl  vs.  a,  Cl  vs.  C^,  and  Cl  vs.  Cq 
variations  for  the  C-141B/AC2  wing  +  fuselage  +  gearpod  configuration. 
Comparing  Figure  29a  wi th  Figure  27a,  we  see  that  the  effect  of  the 
fuselage  and  gearpod  is  to  reduce  the  lift  somewhat  (remember  that  the  fuse¬ 
lage  was  non-metric).  The  reduction  in  lift  is  slightly  higher  at  higher 
angles  of  attack.  The  pitching  moment  data  (Figure  29b)  is  generally  non¬ 
linear,  except  at  lower  values  of  lift,  where  for  a  given  Mach  number,  there 
is  little  variation  in  CM.  The  drag  polars  of  Figure  29c  exhibit  the  same 
sort  of  variation  as  the  wing  alone  configuration,  except  for  somewhat  lower 
lifts  due  to  the  effects  of  the  fuselage. 

Figures  26a  through  26c  summarize  the  force  data  obtained  for  the  C- 1 4 1 B 
baseline  wing  +  fuselage  +  gearpod  configuration.  The  lift  and  pitching 
moment  data  do  not  exhibit  any  unexpected  variations.  The  drag  polars  from 
Figures  26c  and  29c  will  be  used  as  the  basis  for  evaluating  the  C-141B/AC2 
against  the  baseline  C— 1 4 1 . 

Drag  Rise  Characteristics 

Figure  35  gives  the  drag  rise  characteristics  for  the  baseline  C- 1 41 B 
and  the  C - 1 4 I B/AC2  wing  +  fuselage  +  gearpod  configurations.  The  drag 
diverqence  Mach  number  curves  which  were  derived  from  Figure  35,  are  given 
in  Figure  36.  The  original  C- 1 4 1 B/AC2  design  point  is  indicated  in  Figure  36. 
The  target  drag  divergence  Mach  number  at  Cl  *  .6  was  not  achieved  in  the 
wind  tunnel  test,  mainly  because  of  aeroelastic  deformation  of  the  model 
wing,  and  viscous  uncambering  in  the  cove  region.  However,  it  will  be  shown 
later  that  a  very  satisfactory,  high  performance  configuration  was  obtained 
at  CL  ■  .56  and  Mq  -  .78.  * 
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b.  Correlations 


Correlations  of  the  computer  codes  with  the  wind  tunnel  results  were 
made  for  the  C-I41B/AC2  wing  alone,  and  wing  +  pylon-nacelle  configurations. 
The  correlations  were  done  by  trying  to  match  the  measured  wind  tunnel  lift 
by  adjusting  the  angle  of  attack  used  in  the  computer  codes.  The  corre¬ 
lations  were  made  at  nominal  Mach  numbers  of  .78  and  .80.  No  Mach  number 
corrections  were  attempted  because  it  would  have  required  many  more  com¬ 
puter  runs  to  match  lift,  and  because  the  uncertainties  in  predicting  aero- 
elastic  model  deformation  would  overshadow  any  small  Mach  number  correction. 

The  effects  of  model  deformation  under  aerodynamic  loading  can  be  seen 
in  Figure  37,  where  measured  lift  is  compared  to  computed  lift.  The  com¬ 
puted  lift  was  obtained  from  a  full  potential  equation  code  with  iterated 
boundary  layer  (FL022NM) .  Notice  that  when  measured  ordinates  along  with 
twist  values  which  have  been  corrected  for  aerodynamic  loading  are  used,  the 
predicted  Cl  vs.  a  is  much  closer  to  the  data. 

Full  Potential  Equation  Correlations  —  Wing  Alone 

The  code  correlations  shown  in  Figures  38  through  **3  were  made  for  the 
C-14IB/AC2  wing  alone  configuration,  using  a  full  potential  equation  code 
with  iterated  boundary  layer.  The  angles  of  attack  used  in  the  computations 
were  chosen  so  as  to  match  experimental  lift  values  obtained  from  the  force 
balance.  The  correlations  were  done  for  two  Mach  numbers  (-78,  .80)  and 
three  values  of  lift  (CL  *  .33,  -**3,  .60).  For  each  of  these  combinations 
of  Mach  number  and  lift,  we  made  two  calculations.  The  first  calculation 
was  made  using  the  final  design  ordinates,  and  the  design  twist  distribution. 
The  second  calculation  used  measured  model  ordinates  and  estimated  twist. 

After  fabrication,  the  model  was  measured  by  a  special  measuring  machine, 
to  obtain  the  "as  built"  geometry.  Initial  correlations  using  measured 
ordinates  and  design  twist  showed  that  the  twist  distribution  was  wrong. 

The  reason  for  this  was  that  the  model,  having  an  aspect  ratio  of  12,  was 
long  and  thin.  This  tended  to  make  the  model  susceptible  to  aeroelastic 
deformation.  An  additional  factor  was  the  routing  of  140  pressure  tubes. 
Routing  these  tubes  necessitated  cutting  a  large  channel  in  the  lower  surface 
of  the  wing  which  further  weakened  the  structure. 

To  estimate  what  the  aeroelastic  deformation  might  actually  be  when  the 
wing  is  in  the  tunnel  and  under  load,  we  applied  known  loads  to  the  wing  to 
obtain  its  structural  influence  coefficients.  These  influence  coefficients 
together  with  measured  tunnel  pressures  were  used  to  estimate  corrections  to 
the  design  twist  distribution.  The  measured  ordinates  and  estimated  twists 
were  then  used  for  a  second  set  of  correlations-. 

Examining  Figures  38  through  43,  we  see  that  by  using  measured  ordinates 
and  a  twist  distribution  adjusted  for  aeroelastic  deformation,  we  have 
improved  the  correlations.  In  general,  the  best  correlations  are  obtained 
at  the  n  ■  .403  station.  Most  of  the  correlations  show  discrepancies  in  the 
leading  edge  upper  surface,  and  in  the  lower  surface  cove  region.  The  dis¬ 
crepancy  in  the  cove  pressures  is  character i st i c  of  flow  separation  (which 
was  not  modeled) . 
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There  are  several  possible  reasons  for  the  discrepancy  near  the  leading 
edge,  and  we  were  not  able  to  isolate  any  one  reason  in  particular.  The 
leading  edge  discrepancy  is  probably  due  to  several  factors  acting  together. 

The  possibilities  are:  (1)  Correlations  were  made  to  match  total  lift  and 
not  leading  edge  pressures,  (2)  Existence  of  a  short  separated  region  near 
the  leading  edge,  (3)  Transition  strips  which  affect  the  readings  of  some 
of  the  leading  edge  pressure  taps,  (4)  An  actual  twist  distribution  due  to 
aerodynamic  loading  that  is  different  from  the  predicted  twist  distribution, 
and  (5)  The  small  physical  scale  of  the  model  which  is  conducive  to  leading 
edge  irregularities. 

Of  particular  note  in  the  pressure  distributions  is  the  presence  of  a 
shock  wave  with  an  approximately  constant  angle.  This  weak  swept  shock 
wave  was  specified  in  the  target  pressures  used  to  design  the  wing. 

Extended  Smal l  Disturbance  Equation  Correlations  -  -  Wing  Alone 

The  code  correlations  shown  in  Figures  44  and  45  were  made  for  the  C- 1 4 1 B/ 
AC 2  wing  alone  configuration  using  the  extended  small  disturbance  code  ATWP. 
Boundary  layer  effects  were  included  by  using  viscous  linking  process  of  the 
design  procedure.  Angles  of  attack  used  in  the  computations  were  chosen  so 
as  to  match  the  experimental  lift  values  of  the  force  balance.  Correlations 
were  done  for  one  Mach  number  (.80)  and  two  values  of  lift  (Cl=.32,  .49). 

The  small  disturbance  computations  used  design  ordinates  and  twist. 

Examing  Figures  44  and  45  we  see  that  discrepancies  in  the  correlations 
occur  in  the  vicinity  of  the  upper  surface  leading  edge  and  in  the  lower  sur¬ 
face  cove  region.  The  reasons  for  these  discrepancies  are  as  discussed  in 
the  section  on  the  full  potential  equation  correlations. 

Correlations  for  Wing-Pylon-Naoelle  Configuration 

Figure  46  shows  computations  made  using  the  Lockheed  small  disturbance 
wing-pylon-nacelle  code.  Because  of  restrictions  on  the  way  the  nacelle 
computational  grid  is  distributed,  the  computed  and  experimental  results 
shown  in  Figure  46  are  not  exactly  at  the  same  spanwise  location.  But  they 
are  close  enough  to  make  a  meaningful  comparison,  and  thus,  only  the  pressure 
row  n  station  is  given  in  the  figure. 

The  pylon/nacelle  was  positioned  at  n  *  .309*  Although  the  small  dis¬ 
turbance  code  mispredicts  the  wing  leading  edge  upper  surface,  the  pylon/ 
nacelle  interference  predictions  are  excellent.  At  locations  inboard  of  the 
nacelle  the  interference  effects  are  strong.  Outboard  of  the  nacelle  the 
interference  is  minor.  The  predictions  follow  the  experimental  trends 
closely. 

c.  Effect  of  Wind  Tunnel  Walls 

Figure  47  shows  the  effect  on  computed  wing  pressures  of  including 
measured  wind  tunnel  wall  pressures  as  boundary  conditions.  Two  computed 
sets  of  pressures  are  shown.  The  first  is  a  small  disturbance  solution  at 
M  ■  .8  and  angle  of  attack  »  .7°  with  free  air  boundary  conditions.  The 
second  is  a  small  disturbance  solution  where  measured  wind  tunnel  wall 
static  pressures  were  imposed  as  boundary  conditions.  The  biggest  effect 
of  the  wind  tunnel  walls  is  to  move  the  upper  surface  shock  forward.  The 
effects  on  the  lower  surface  pressures  are  not  large,  except  at  the  two 
inboard  stations. 
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d.  Evaluation  of  Pressure  and  Force  Balance  Data 

We  took  the  measured  pressures  on  the  C— 1 4 1 B/AC2  model,  integrated  them 
to  obtain  section  lifts,  and  then  integrated  the  section  lifts  to  obtain 
total  lift.  We  then  compared  the  total  lift  obtained  from  the  pressure 
integration  with  lift  obtained  from  the  force  balance.  In  general,  the  force 
balance  lift  was  higher  than  the  pressure  lift.  A  difference  between  the 
two  lifts  is  to  be  generally  expected  and  would  normally  not  exceed  5  or  6 
percent.  In  our  case,  however,  differences  of  the  order  of  10  or  12  percent 
were  obtained  for  some  runs.  Figure  A8  illustrates  the  problem. 

The  question  then  is  what  would  account  for  such  a  discrepancy.  Either 
the  force  balance  readings  were  incorrect,  or  there  was  something  wrong  with 
the  results  from  the  pressure  integration.  We  checked  the  force  balance 
electronics  and  data  acquisition  system,  and  found  them  to  be  operating 
properly.  As  a  further  check  on  the  force  balance  we  installed  a  Lockheed 
ATA-V  wing  model  in  the  wind  tunnel,  tested  it,  and  compared  these  results 
with  those  obtained  from  a  previous  tunnel  entry;  the  lift  and  drag  data  from 
the  two  entries  were  consistent.  We  then  integrated  ATA-V  surface  pressures 
and  compared  them  to  the  force  data.  These  comparisons  showed  differences 
that  were  in  the  acceptable  range.  It  can,  therefore,  be  concluded  that 
neither  the  force  balance  readings  nor  the  pressure  integration  methodology 
are  responsible  for  the  observed  C  —  1 4 1 B/AC2  lift  discrepancies. 

We  then  carefully  checked  the  pressure  data,  eliminating  or  fairing 
through  suspicious  pressure  readings,  such  as  those  in  the  vicinity  of  the 
transition  strips.  By  integrating  these  new  curves,  we  could  eliminate  about 
half  of  the  excessive  discrepancy.  Although  we  could  not  pinpoint  the  exact 
cause  of  the  remaining  portion  of  the  discrepancy,  we  have  identified  a 
number  of  possible  contributing  factors.  These  include  a)  the  existence 
of  a  small  region  of  separated  flow  on  the  upper  surface  leading  edge,  b) 
the  aeroelastic  twist  and  deflection  due  to  the  flexible  wing  structure,  and 
c)  the  viscous  lag  in  the  smal ler- than -normal  (.020")  tubing  used  because  of 
the  relatively  small  size  of  the  wing.  The  latter  (lag)  effects  were  re¬ 
checked  experimentally,  and  were  found  to  account  for  a  pressure  error  of 
less  than  one  percent  under  even  the  most  adverse  sampling  conditions. 

In  summary,  then,  our  investigations  indicate  that  the  force  balance 
readings  of  lift  are  correct.  Because  of  a  combination  of  factors  which 
adversely  affect  the  precision  of  our  surface  pressure  readings,  however, 
pressure-integrated  lift  values  cannot  be  expected  to  produce  an  exact  match 
for  comparable  force  measurements. 
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5.  TASK  5  '  EVALUATION  OF  DESIGNS 


There  were  two  subtasks  associated  with  Task  V.  The  first  subtask  was 
to  use  the  results  of  the  wind  tunnel  tests  to  quantify  the  performance 
improvements  of  the  new  C-141B/AC2  configuration  over  the  baseline  C- 1 4 1 B , 
and  in  addition,  to  determine  the  cost  savings  attributable  to  the  new 
design  procedure.  The  second  subtask  was  to  evaluate  the  design  procedure. 


Performance  Improvements 


The  C-1418/AC2  design  conditions  and  wing  geometry  are  substantially 
different  from  the  C- 1 4 ! .  Consequently,  a  comparison  of  wing  aerodynamics 
does  not  by  itself  provide  a  true  measure  of  the  effectiveness  of  the  new 
wing.  For  example,  a  comparison  of  drag  polars  at  the  C- 1 4 1 B/AC2  design 
Mach  number  would  not  be  meaningful  because  the  C— 1 4 1  was  designed  to  cruise 
at  .77  Mach,  and  its  wing  is  well  into  drag  rise  at  .80  Mach.  Consequently, 
the  efficiency  of  the  design  procedure  is  evaluated  by  comparison  of  complete 
airplane  performance  capabilities  rather  than  by  reference  to  incremental 
aerodynamic  characteristics.  To  make  these  comparisons,  flight  aerodynamics 
for  the  C- 1 4 1 B/AC2  are  extrapolated  from  wind  tunnel  data  using  the  known 
C - 1 4 1  flight  characteristics  as  a  calibration. 


In  the  Compressible  Flow  Wind  Tunnel,  the  tunnel  top  and  bottom  walls 
are  relatively  close  -  approximately  3i  mean  chords  -  to  the  wing.  Although 
a  procedure  for  taking  wall  effects  into  account  is  under  development,  it 
is  not  at  the  point  where  it  can  De  used  with  confidence.  Accordingly,  the 
analysis  method  we  used  is  based  on  a  comparison  of  uncorrected  measured 
and  estimated  drag  differences  for  the  two  wing-fuselage-gear  pod  config¬ 
urations.  The  drag  estimation  technique  is  known  to  agree  well  with  C- 1 4 1 
f 1 ight  experience. 


A  drag  estimation  at  M  *  0.70  was  made  for  each  configuration  at  the 
wind  tunnel  Reynolds  number  and  then  compared  with  the  measured  model  drag. 
The  resulting  drag  differences  are  shown  in  Figure  49.  The  variation  with 
lift  coefficient  is  identical  for  the  two  designs,  and  the  C-I41B/AC2  drag 
is  10  counts  less  than  the  C-I4IB  drag.  The  actual  drag  of  the  C-141B  full" 
scale  aircraft  is  known  from  flight  tests  and  is  reproduced  by  the  drag 
estimation  method  employed.  Assuming  that  the  drag  increment  is  insensitive 
to  Reynolds  numbers,  the  full-scale  drag  of  the  C-141B/AC2  aircraft  at 
M  ■  0.70  can  be  obtained  by  subtracting  10  counts  from  the  full-scale  esti¬ 
mation  drag  polar  for  the  C— 1 4 1 B . 


The  next  step  in  the  drag  analysis  procedure  was  the  determination 
of  the  drag  rise  Mach  number  and  the  compressibility  drag  increment.  A 
direct  comparison  of  the  measured  drag  rise  characteristics  of  the  wing- 
fuselage-gear  pod  configurations  at  constant  Cj_  is  presented  in  Figure  35. 
The  corresponding  drag  divergence  Mach  numbers,  Mg,  are  shown  as  a  function 
of  lift  coefficient  in  Figure  36.  The  target  Mg  value  for  the  C-141B/AC2 
design  of  M  *  0.80  at  Cj_  ■  0.60  was  not  achieved.  This  results  partly  from 
the  viscous  uncambering  of  the  cove  region  at  the  relatively  low  tunnel 
Reynolds  number,  and  partly  from  the  aeroelastic  deformation  of  the  model 
wing  at  the  high  dynamic  pressures  of  the  Compressible  Flow  Wind  Tunnel. 


Reducing  the  full-scale  design  Mach  number  and  lift  coefficient  to  M  ■  O.78 
and  ■  0.56  in  accordance  with  the  measured  data  of  Figure  36,  nevertheless 
results  in  a  highly  satisfactory  advanced  technology  C-141B/AC2  aircraft 
design.  Figure  50  gives  comparisons  of  payload-range  performance.  These 
comparisons  show  the  actual  payload-range  performance  of  the  C-141B,  the 
C-141B/AC2  target  performance  which  was  based  on  what  performance  improve¬ 
ments  we  would  expect  to  obtain  by  applying  the  new  design  method,  and  the 
predicted  (based  on  wind  tunnel  test  results)  C-141B/AC2  performance  which  we 
actually  did  obtain.  Of  note  is  the  better  than  targeted  ferry  range  of  the 
C-141B/AC2  made  possible  by  the  thick  wing  resulting  from  the  application  of 
the  new  wing  design  method. 

The  reduction  in  aircraft  weights  and  fuel  made  possible  by  advanced 
technology  are  shown  in  Figure  51.  The  predicted  gross  weight  reduction  is 
21%.  This  is  1%  better  than  the  target  value  of  20%.  The  predicted  empty 
weight  reduction  is  24%,  which  is  4%  better  than  the  target  value  of  20%. 

The  predicted  block  fuel  decrease  is  36%  which  is  equal  to  the  target  value 
of  36%.  Thus,  with  the  exception  of  the  .80  cruise  Mach  number,  the  study 
design  objectives  have  been  achieved. 

b.  ATT  Cost  Benefit  Analysis 

The  question  is  always  asked:  "Can  you  quantify  the  benefits  of  new 
computational  methods"?  In  this  section  we  will  do  that  for  the  ATT  design 
procedure. 

Previous  Method 

During  1977  and  the  early  part  of  1978,  a  transonic  wing  was  designed 
at  Lockheed-Georgia  for  the  Advanced  Transport  Aircraft  (ATA)  configuration. 
The  ATA  configuration  was  a  company- funded  design  exercise  to  oroduce  an 
Advanced  Transport  Cargo  Aircraft.  The  wing  was  designated  ATA  -  V. 

The  design  process  for  the  ATA  -  V  wing  involved  a  total  of  7  different 
computer  programs  and  7  engineers.  The  computer  programs  used  were: 

(1)  Carlson  Transonic  inverse  2-D  code 

(2)  TAP  -  2-D  analysis  code 

(3)  WPRESS  -  2-D  analysis  code 

(4)  WPRESS2  -  improved  2-D  analysis  code 

(5)  L-7  3~0  lifting  surface  analysis  code 

(6)  ATWP  -  3_D  Transonic  analysis  code 

(7)  AFSM00TH  -  interpolation  code 

COMPUTER  COSTS 


The  following  runs  were  made  for  the  ATA  -  V  Design: 


(1) 

CODE 

CARLSON 

NO.  OF  RUNS 

20 

COST 
$  600 

(2) 

TAP 

35 

875 

(3) 

WPRESS  &  WPRESS2 

30 

150 

(4) 

L*7 

20 

300 

21 


(5) 

ATWP 

(6) 

AFSM00TH 

24  6,000 

Charges  not  significant 


TOTAL  DOLLARS 

$7,925 

LABOR 

COST 

2,196  Manhours 

$73,675 

TOTAL  ATA-V  DESIGN 

COST 

$81 ,600 

Present  Method 

The  numerical  optimization  design  procedure  developed  under  the  ATT 
contract  was  used  in  the  high  speed  wing  design  of  the  projected  new 
military  transport  currently  known  as  the  C-X.  This  gives  us  a  cost  figure 
for  a  wing  design  exclusive  of  any  design  procedure  development  costs,  which 
is  what  we  need.  The  C-X  wing  design  involved  3  engineers  over  a  span  of 
16  days.  Here's  the  cost  breakdown: 


COMPUTER  RUNS  COST 

32  $35,000 

LABOR  COST 

420  manhours  $14,091 

TOTAL  COSTS:  $49,091 


In  looking  at  the  relative  costs,  we  can  see  that  the  ATT  procedure 
required  $27,075  more  dollars  in  computer  charges,  but  $49,584  less  in 
labor,  thus  resulting  in  a  net  savings  of  $32,509.  This  amounts  to  a  hot 
reduction  of  total  wing  design  costs. 

Such  savings  are  significant  in  terms  of  dollars.  What  is  even  more 
significant  is  that  the  calendar  time  required  to  do  the  design  is  substan- 
ially  less.  The  ATA-V  design  was  spread  over  a  period  of  14  months.  If  a 
concentrated  effort  was  made,  the  ATA-V  design  couid  have  been  done  in  six 
months  under  the  old  procedure.  This  means  that  using  the  new  design  pro¬ 
cedure  would,  in  a  case  like  this,  decrease  transonic  wing  design  from  six 
months  to  just  over  two  weeks. 
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The  potential  for  even  further  reduction  in  design  cost  is  inherent 
in  the  ATT  design  procedure.  Under  Lockheed-Georg ia  internal  funding  we  have 
carried  out  proof  of  concept  studies  where  2-D  airfoil  optimization  computer 
time  has  been  decreased  by  a  factor  of  5.  For  3-0  wing  design  using  num¬ 
erical  optimization  the  factor  of  5  is  probably  a  bit  conservat ive.  Never¬ 
theless,  with  improvement  by  a  factor  of  5  in  a  design  such  as  the  C-X  the 
computer  charge  would  then  be  $7,000  for  a  total  including  labor  of  $21,091. 
This  is  26%  of  the  ATA-V  design  cost. 

c.  Critique  of  Design  Procedure 


We  have  shown  that  using  numerical  optimization  is  a  successful  way  to 
design  wings  which  produce  desired  cruise  pressure  distributions  and  that 
the  method  can  be  incorporated  within  the  framework  of  current  aircraft 
design  procedures.  However,  we  have  identified  several  areas  where  the 
current  design  procedure  may  be  improved.  These  areas  are: 

1 .  Design  variables 

2.  User  expertise 

3.  Computation  time 


Design  Variables 

Sine  deformation  shape  functions  of  the  form  sinn(irx)  provide  the 
primary  means  of  modifying  the  wing  geometry  to  produce  the  desired  pressures 
The  intent  of  these  functions  i^  to  provide  local  geometric  control.  For 
example,  the  shape  function  sm  ?rx  oroduces  maximum  geometric  change  at  the 
50%  chord  station  as  shown  in  Figure  52..  Also  shown  in  that  figure  is  the 
change  in  curvature  produced  by  the  shape  function.  Although  the  maximum 
curvature  change  occurs  at  50%  chord,  there  are  two  other  locations  of  signif 
icant  curvature  change.  Since  curvature  plays  a  dominant  role  in  the  devel¬ 
opment  of  the  flow  field,  the  non-local ized  curvature  changes  produced  by 
the  sine  shape  functions  can  cause  undesired  changes  in  the  flow  field  and 
that  introduces  ambiguity  in  the  optimization  process. 


A  remedy  for  this  is  to  use  design  variables  based  on  the  second 
derivative  of  the  airfoil  surface.  Such  shape  functions  not  only  localize 
the  curvature  variations,  but  also  produce  very  smooth  airfoils.  Research 
is  presently  underway  at  Lockheed-Georg ia  and  at  the  Ames  Research  Center, 
NASA,  to  develop  efficient  curvature-based  design  variables. 


User  Expertise 

Engineers  who  are  skilled  aerodynamici sts  and  who  are  familiar  with 
numerical  optimization  are  needed  to  successfully  use  optimization  in 
aircraft  design.  This  need  exists  because  of  0)  the  requirement  to 
accurately  specify  desirable  pressure  distributions  which  produce  realistic 
wing  geometries,  and  (2)  difficulties  encountered  in  selecting  design 
shape  functions  which  will  produce  the  desired  flow  field  modifications.  The 
latter  difficulty  can  be  ameliorated  by  the  use  of  curvature-based  design 
variables. 
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Clearly,  the  difficulties  encountered  in  the  specification  of  desirable 
pressures  are  accentuated  for  multiple  design  point  aircraft  such  as  super¬ 
sonic  crui se/transon ic  maneuver  fighters.  For  transport  aircraft,  consid¬ 
eration  of  wing  weight  and  drag  reduction  must  be  balanced  against  one 
another.  One  possible  solution  to  this  problem  is  to  conduct  studies 
to  identify  sensible  and  desirable  pressure  distributions  for  different 
missions. 

Another  alternative  is  to  use  design  objectives  based  on  aero¬ 
dynamic  forces  and  moments.  For  this  approach  to  be  successful,  the 
accuracy  of  computed  aerodynamic  forces  and  moments,  in  particular  drag, 
must  be  improved.  Experience  with  current  numerical  aerodynamic  methods, 
even  in  two  dimensions,  has  shown  that  inaccuracies  in  drag  calculations 
can  make  realistic  and  reliable  numerical  optimization  difficult.  If 
sufficiently  consistent  and  accurate  drag  calculation  techniques  can  be 
developed,  then  the  use  of  design  objectives  based  on  integrated  aerodynamic 
parameters  would  best  take  advantage  of  the  capabilities  offered  by  numer¬ 
ical  optimization. 

Computation  Time 

Between  5  to  10  hours  of  computation  time  on  a  CDC  7600-class  computer 
are  needed  to  perform  a  wing  design  using  the  current  numerical  optimization 
scheme.  These  relatively  large  times  are  caused  by  the  multitude  of  non¬ 
linear  aerodynamic  solutions  required  during  the  optimization  process.  The 
three  obvious  ways  of  reducing  the  computation  times  are  (1)  use  more 
efficient  computers,  (2)  use  better  solution  algorithms,  and  (3)  reduce  the 
number  of  non-linear  solutions. 


The  first  two  ways  are  related  and  they  involve  the  use  of  new 
algorithms  such  as  approximate  factorization  schemes  on  new  vector  computers 
such  as  the  CRAY-1  and  the  CDC  Cyber  203.  Significant  research  is  being 
devoted  to  this  task. 

The  third  way  can  be  approached  by  at  least  two  means.  One  approach 
is  to  develop  a  versatile  and  reliable  3-D  inverse  transonic  method  in  which 
the  wing  geometry  is  computed  directly  from  the  specified  pressures.  Such 
a  method  would  require  about  the  same  computation  time  used  in  transonic 
flow  analysis.  A  deficiency  in  this  approach  is  that  constraints  are  diffi¬ 
cult  to  impose.  Nevertheless,  an  inverse  method  could  produce  a  wing  that 
is  nearly  an  acceptable  design.  Numerical  optimization  could  then  be  used 
for  the  final  design  refinements.  The  geometric  changes  might  be  expected 
to  be  less  than  for  a  complete  optimization  design,  and  fewer  design  vari¬ 
ables  might  be  required.  Thus,  the  number  of  non-linear  solutions  needed 
in  the  optimization  process  should  be  reduced. 

The  second  approach  to  reducing  the  number  of  non-linear  solutions  in 
fact  involves  the  replacement  of  fine  grid  solutions  with  coarse  grid 
results.  To  maintain  accuracy,  the  coarse  grid  results  are  corrected  to 
equivalent  fine  grid  accuracy  using  Nixon's  strained  coordinate  scheme  18, 

The  possibility  of  such  an  approach  has  already  been  explored  by  lockheed- 
Georgia  and  Nielsen  Engineering  and  Research  scientists.  In  a  proof  of 
concept  study,  optimization,  together  with  Nixon's  strained  coordinate 
scheme  have  been  coupled  to  a  2-D  transonic  airfoil  code.  With  this  scheme, 
computer  time  for  airfoil  design  using  numerical  optimization  has  been 


reduced  by  close  to  a  factor  of  5- 


By  implementing  a  design  procedure  incorporating  an  inverse  method 
together  with  numerical  optimization  and  strained  coordinates,  we  can  expect 
to  reduce  wing  design  computation  time  from  the  current  5  to  10  hour  range 
to  approximately  1  to  2  hours  on  a  CDC  7600.  By  using  an  algorithm  which 
takes  advantage  of  new  vector  processing  computers,  computer-aided  wing 
transonic  aerodynamic  design  in  less  than  1/2  hour  can  be  forecast. 

6.  TASK  6  -  FINAL  DESIGN  PROCEDURE 

The  wing  design  procedure  involves  the  use  of  combinations  of  the 
following  computer  programs: 

o  Full  Potential  Inviscid  NCR 

o  Extended  Small  Disturbance  NCR  and  FCR 

o  Numerical  Optimization 

o  Two-Dimensional  Boundary  Layer 

o  Linking  Codes 

o  Subsonic  Panel  Code 

o  Transonic  Wing/Pylon/Nacelle  Code 

The  steps  involved  in  the  procedure  are  shown  in  Figure  53. 

a.  Starting  Wing 

Starting  wing  geometry  is  not  critical  to  the  success  of  the  method; 
however,  computation  costs  can  be  reduced  by  selecting  a  wing  with  well- 
behaved  flow  near  the  design  condition.  For  transonic  transports  the  air¬ 
foil  sections  designed  for  this  contract  should  be  satisfactory  for  typical 
wing  planforms.  Alternatively,  a  wing  with  which  the  user  has  some  previous 
experience  can  be  used,  or  a  wing  with  an  airfoil  section  designed  with  2-D 
design  methods  (e.g.,  ref.  12,  17)  might  be  employed. 

The  wing  geometry  is  defined  by  a  number  of  spanwise  geometric  control 
stations  connected  by  spanwise  linear  loft  elements.  We  have  found  that 
this  approach  provides  sufficient  design  flexibility  and  at  the  same  time 
produces  a  wing  which  is  easy  to  manufacture  by  avoiding  compound  curvature. 
Our  experience  shows  that  a  four  control  station  wing  (root,  break,  mid¬ 
semispan,  tip)  is  adequate.  However,  an  additional  station  between  the  root 
and  break  might  be  needed  in  some  applications. 

b.  Inviscid  Solution 


At  this  point,  the  user  should  select  the  transonic  analysis  code  -  FPE 
or  ESD  NCR/FCR.  Based  on  our  experience,  FPE  optimization  is  more  reliable, 
but  ESD  NCR  might  be  satisfactory  if  the  wing  leading  edge  is  not  being 
changed.  Although  we  have  not  done  a  conclusive  analysis,  we  believe  ESD/FCR 
design  should  be  avoided. 
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Having  selected  the  analysis  code,  compute  an  inviscid  solution  for  the 
starting  wing  at  the  design  Mach  number.  The  angle  of  attack  should  be 
selected  to  produce  a  lift  coefficient  which  is  about  10%  greater  ttwi  the 
design  untrimmed  lift  coefficient.  The  geometry  file  must  be  cataloged  for 
later  use  and  the  pressure  distribution  must  be  cataloged  for  use  in  the 
boundary  layer  analysis. 


Attach  the  pressure  distribution  data  file  and  use  the  2-D  boundary 
layer  code  to  compute  displacement  thicknesses  at  the  wing  geometric  control 
stations.  At  this  point,  the  user  should  select  the  design  Reynolds  number. 
If  the  wing  is  being  designed  for  technology  demonstration,  tunnel  Reynolds 
number  should  be  used.  Flight  Reynolds  number  is  advisable  for  aircraft 
design  applications. 

Catalog  the  displacement  thickness  file  for  use  in  the  next  step, 

d.  Link  Up 

Using  the  Link-Up  program  with  the  geometry  and  displacement  thickness 
files,  compute  the  starting  viscous  pressure  distribution. 

NOTE:  The  inviscid  solution,  boundary  layer  analysis,  and  link-up 
steps  could  be  replaced  by  a  3_D  viscous  analysis  (e.g.  the 
viscous  FPE  code  available  from  NASA-Ames).  However,  this 
approach  can  introduce  some  inconsistency  in  the  design  pro¬ 
cedure. 


e.  Design  Pressures 

Design  pressures  are  specified  near  each  geometric  control  station.  The 
tip  pressures  should  be  defined  near  the  90%  span  station  to  avoid  inaccuracies 
in  solutions  near  the  tip.  Our  experience  shows  that  design  pressures  with  a 
single  weak  shock  wave  outboard  on  the  wing  should  be  selected.  The  root 
upper  surface  pressures  should  have  a  mild  re -compress ion  gradient.  For  many 
applications,  the  upper  surface  pressure  distributions  used  in  this  study 
should  be  quite  satisfactory.  The  lower  surface  pressures  should  be  selected 
to  produce  the  desired  aft  loading  and  wing  thickness. 

The  subsonic  panel  code  is  used  to  ascertain  the  fusel age/gearpod  inter¬ 
ference  effects.  These  should  be  accounted  for  in  the  design  pressures.  The 
interference  effects  of  the  pylon/nacelle  installation  can  be  determined  by 
running  the  transonic  wing-pylon/nacelle  code. 

f .  Des i gn 

Wing  design  is  accomplished  by  designing  the  upper  surface  first  and 
then  the  lower  surface.  Upper  surface  design  is  done  station  by  station 
from  root  to  tip.  Specifically,  the  root  station  is  designed  to  produce 
specified  root  pressures  while  keeping  other  stations  fixed.  The  centerline 
station  wing  carry-through  geometry  is  automatically  changed  to  be  the  same 
as  the  root  station.  We  have  found  that  it  is  best  to  change  the  entire  airfoil 
upper  surface  at  a  control  station  if  the  pressures  at  that  control  station  are 
shock  free.  Otherwise,  the  design  variables  which  affect  geometric  changes 
upstream  of  the  shock  should  be  used  to  produce  the  desired  pressures  ahead 
of  the  shock.  Next,  the  recompression  downstream  of  the  shock  is  designed. 
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Once  the  pressure  distribution  at  the  root  is  satisfactory,  the  root  air¬ 
foil  geometry  is  fixed  at  the  geometry  that  came  out  of  the  root  design, 
and  the  next  outboard  station  is  designed. 

After  each  upper  surface  station  is  designed,  the  lower  surface  is 
designed  proceeding  from  root  to  tip.  The  lower  surface  pressures  should 
be  selected  to  produce  the  desired  wing  thickness  and  aft  loading.  The 
entire  lower  surface  at  a  span  station  can  be  designed  simultaneously. 

After  the  entire  wing  has  been  modified,  the  procedure  can  be  repeated, 
as  required.  Our  experience  shows  that  two  iterations  of  the  wing  design 
procedure  are  generally  required. 

g.  Invi sci d  Anal ys i s 

The  designed  fluid  wing  is  next  analyzed  using  the  inviscid  analysis 
code  which  was  used  during  design.  The  resulting  pressures  and  the  wing 
geometry  are  cataloged  for  use  in  the  next  two  steps. 

h .  Boundary  Layer  Analysis 


The  2-D  boundary  layer  code  is  used  to  compute  displacement  thickness 
distributions  at  each  control  station.  The  output  of  the  analysis  should 
be  carefully  reviewed  to  define  separation  locations.  Separation  in  the 
cove  should  be  avoided,  but  separation  near  the  trailing  edge  (~37%  chord) 
might  be  unavoidable  and  is  not  a  serious  problem.  The  user  should  use  all 
separation  indicators  (d6*/dx,Cf,  form  factor)  in  making  an  engineering 
judgment  to  define  separation. 

The  displacement  thickness  file  is  cataloged  for  use  in  the  next  step, 

i.  Link-Down 


The  link-down  program  subtracts  the  5*  from  the  fluid  geometry  to  produce 
the  solid  wing.  The  user  should  carefully  review  the  resulting  solid  wing 
geometry  for  smoothness.  We  have  found  that  one  or  two  smoothings  are  desirable. 

j .  Viscous  Analysis 

The  last  step  of  the  procedure  is  a  viscous  analysis  of  the  solid  wing. 

The  analysis  should  be  done  using  both  the  linking  technique  of  the  design 
method  and  a  coupled  viscous  analysis  method.  The  results  of  this  analysis 
should  compare  closely  to  the  results  of  the  inviscid  analysis  of  the  fluid 
wi  ng. 

k.  A  Final  Note 


The  design  procedure  involves  a  great  deal  of  user  interactions  with  the 
computer  codes.  This  was  done  on  purpose.  Engineering  judgment  is  required 
in  many  key  steps,  and  in  particular,  steps  involving  boundary  layer  analyses. 
We  do  not  believe  that  our  ability  to  compute  separation,  transition,  and 
equivalent  fluid  geometry  is  at  present,  sufficiently  reliable  to  permit 
hands-off  design. 
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SECTION  I  I  I 


CONCLUDING  COMMENTS 

We  have  developed  a  new  transonic  wing  design  method  using  the  numerical 
optimization  scheme.  We  have  also  shown  that  new  computational  methods  offer 
a  means  for  the  aerodynamic  design  of  wings  with  transonic  performance  super¬ 
ior  to  that  which  could  be  obtained  using  previous  design  techniques.  The 
method  is  relatively  easy  to  use,  and  it  is  compatible  with  established 
industry  design  procedures.  By  using  the  new  method,  a  40%  to  50%  reduction 
in  the  cost  associated  with  wing  cruise  aerodynamic  design  is  obtainable. 

By  incorporating  the  strained  coordinate  scheme  into  the  method,  cost  re¬ 
ductions  should  approach  75%.  Additionally,  adapting  a  strained  coordinate 
version  of  the  method  to  the  new  vector  processing  computers  should  make 
possible  efficient  wing  design  using  less  than  i  hour  computation  time. 


REFERENCES 


1.  Boppe,  C.  W. :  "Calculation  of  Transonic  Wing  Flows  by  Grid  Embedding," 

AIAA  Paper  No.  77_207,  January  1977. 

2.  Boppe,  C.  W. :  "Computational  Transonic  Flow  About  Realistic  Aircraft 
Configurations,"  AIAA  Paper  No.  78-104,  January  1978. 

3.  Mason,  W.  H. ,  et  al.:  "An  Automated  Procedure  for  Computing  the  Three- 
Dimensional  Transonic  Flow  Over  Wing-Body  Combinations,  Including 
Viscous  Effects,"  AFFDL-TR-77-122,  October  1977. 

4.  Jameson,  A.:  "Iterative  Solution  of  Transonic  Flows  Over  Airfoils  and 
Wings,"  Corrn.  Pure  Appl.  Math.,  Volume  27,  1974,  pp.  283-309. 

5.  Jameson,  A.;  and  Caugbey,  D.  A.:  "A  Finite  Volume  Method  for  Transonic 
Potential  Flow  Calculation,"  AIAA  Paper  No.  77-635,  June  1 977 - 

6.  Blackerby,  W.  T. ;  and  Smith,  P.  R.:  "Aerodynamic  Investigation  of  C— 1 4 1 
Leading  Edge  Modification  for  Cruise  Drag  Reduction,"  AFFDL  TR79-3059, 

June  1979. 

7.  Caughey ,  D.  A.;  and  Jameson,  A.:  "Accelerated  Iterative  Calculation 
of  Transonic  Nacelle  Flow  Fields,"  AIAA  Journal,  Volume  15,  No.  10, 

October  1977,  PP<  1474-1480. 

8.  Vanderplaats ,  G.  N. ;  "CONMIN  -  A  Fortran  Program  for  Constrained 
Function  Minimization,"  NASA  TM  X-62282,  August  1973. 

9.  McNally,  W.  D. ;  "FORTRAN  Program  for  Calculating  Compressible  Laminar 
and  Turbu  lent  Boundary  Layers  in  Arbitrary  Pressure  Gradients,"  NASA 
TN  0-5681,  1970. 

10.  Bartlett,  D.  W. :  "Wind  Tunnel  Investigation  of  Several  High  Aspect  Ratio 
Supercritical  Wing  Configurations  on  a  Wide-Body-Type  Fuselage,"  NASA  TM 
X-71996,  1977. 

11.  Kline,  J.:  "Wind  Tunnel  Tests  of  a  0.030  Scale  Full-Span  Model  of  the 
Lockheed-Georgia  ATA  Series  5  Transport,"  Calspan  Rpt.  AA-4007-W-7,  1978. 

12.  Carlson,  L.  A.:  "Transonic  Airfoil  Analysis  and  Design  Using  Cartesian 
Coordinates,"  AIAA  Proceedings,  Second  Computational  Fluid  Dynamics 
Conference,  Hartford,  Conn.,  1975. 

13.  Bauer,  F.,  et.  al.,  "Supercritical  Wing  Section,"  Lecture  Notes  in 
Economics  and  Mathematical  Systems  66,  Spri nger-Verlag,  1972. 

14.  Jareson,  A.:  "Iterative  Solution  of  Transonic  Flows  Over  Airfoils  and 
Wi.igs,"  Comm.  Pure  and  Appl.  Math,  Volume  27,  1974,  pp.  283-309. 

15.  Pounds,  G.  A.;  and  Stanewsky,  E.:  The  Research  Compressible  Flow  Facility. 
Lockheed-Georgia  Company  ER-9219,  1967. 

16.  Braslow,  Albert  L. ;  and  Knox,  Eugene  C.:  Simplified  Methods  for  Determining 
of  Critical  Height  of  Distributed  Roughness  Particles  for  Boundary  Layer 

29 


m 


TABLE  2:  C-141  PRESSURE  TAP  LOCATIONS 
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TABLE  3:  C-I41/AC2  PRESSURE  TAP  LOCATION 
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TABLE  **:  C-I41B/AC2  MEASURED  MODEL  ORDINATES  6  ESTIMATED  TWISTS  (CONTINUED) 


CFVT052  -  C-141B  CONFIGURATION 


TABLE  6:  RUN  SUMMARY  TEST  CFWT05? 
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TABLE  6:  RUN  SUMMARY  TEST  CFWT053  C-I41B/AC2  CONFIGURATION  (CONTINUED) 
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TABLE  6:  RUN  SUMMARY  TEST  CFWT053  C-141B/AC2  CONFIGURATION  (CONTINUED) 


TABLE  7:  RUN  SUMMARY  TEST  CFWT052 
C-HtlB  CONFIGURATION 


CONTRACT  RESEARCH 


PRODUCTS 


ON-GOING  COMPANY  RESEARCH 


C-141D/AC2  SELECTED  TRANSPORT  COLOURATION 
RANGE:  3500.  N.M.  PATLOAUx  73000.  CRUISE  MACIIs  .000  DELTA  MACH:  .000 

CENTER  SECTION  FUEL  MAS  USED  IN  Y-RANGE  CALCULATION.  HI5C.  ORAGs  17.0  COUNTS. 
.0  PERCENT  2-NO  SEGMENT  L/D 

THRUST  ETA  H— 1000  M/S  CL  AR  SWEEP  T/C  SW  W-RAMP  OW  CR-L/D  FVR 
38774.  .80  33.  130.00  .36312.00  23.00  .1087  1960.  261118.  118273.  21.17  1.039 

/—DESIGN  POINT  DATA—/ 

W-RAMP  OOC/SM  WTTO  WTLNO  FAA1L  LOIST  T/W  ORAO  CLMXT  CLMXL  FLPT  FLRL 

261118.  .000  261118.  227196.  7512.  4113.  .2970  .0423  1.76  2.23  10.0  40.0 
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Figure  5.  c- 1 4 1 B/AC2  Configuration  Data  and  Target  Performance 
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GEOMCTftIC  SUMMARY 


Figure  5-  Continued 
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Figure  5.  Concluded 
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RANGE  -  1000  NM 


Figure  6.  Payload/Range  Comparison 
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Figure  7.  Impact  of  Advanced  Technologies  on  the  C- 1 ^ 1 B/AC2 
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Figure  10.  Transport  Wing  Design  Procedure 
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Figure  15.  Continued 
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Figure  15.  Concluded 


Figure  16.  Continued 
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Figure  16.  Concluded 


Figure  22.  Lockheed-Georgia  Compressible  Flow  Wind-Tunnel 


23.  Photograph  of  Test  Apparatus 
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Figure  27a.  vs.  a  C-141B/AC2  Wing  Alone 
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Figure  23c.  C(_  vs.  Cq  C-141/AC2  Wing  +  Pylon/Nacelle 
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Figure  30a.  Cj_  vs.  •  C-141B/AC2  Wing  +  Fuselage  +  Gearpod  +  Pylon/Nacei le 
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iFigure  32a.  Cj_  vs.  a  C-141B/AC2  Wing  +  Fuselage 
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Figure  3^.  C-I^IB/AC  Wing  Alone  Isometric  Pressure  Plots 
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Figure  3^*  (Continued)  (a  -  0°) 


112 


CONDITIONS 


MACH  NO 

-  0.802 

ALPHA 

0.483 

CL 

0.505 

CD 

0.025 

CM 

0.015 

CLG 

0.505 

Figure  34.  (Continued)  (a  -  5°) 


Figure  34.  (Continued)  (a  -  1.5°) 
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Figure  31*.  (Concluded)  (a  -  3°) 
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Figure  35.  Measured  Drag  Rise  Characteristics 
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Figure  36.  Measured  Drag  Divergence  Mach  Numbers 
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Figure  38.  (Continued)  (4) 


Figure  38.  (Concluded)  (5) 
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Figure  39.  Full  Potential  Correlation:  M  -  0.78,  CL  - 
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Figure  kO.  (Continued)  (2) 
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Figure  40.  (Concluded)  (5) 
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Figure  42.  (Continued)  (3) 


CFVT053  RUN  35 


FL022NM  DESIGN  ORDINATES 
AND  TWIST 

FL022NM  MEASURED  ORDINATES 
AND  PREDICTED  TWIST 


n  =  .653 


o  o  o 


CFVT053  RUN  35 


— FL022NM  DESIGN  ORDINATES 
AND  TWIST 

-FL022NM  MEASURED  ORDINATES 
AND  PREDICTED  TWIST 


ri  =  .  861 


Figure  43-  (Concluded)  (5) 
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Figure  44.  Continued  (2) 
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Figure  44.  Continued  (3) 
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Figure  1*7.  Effect  of  Wind  Tunnel  Wall  Boundary  Conditions, 
M  -  0.80,  o  -  0.7°. 
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Figure  48.  Pressure  Integration  vs.  Force  Balance  Comparisons 
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Figure  49.  Comparison  of  Measured  Model  Drag  With  Estimates 
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Figure  53 •  Design  Procedure 
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